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ABSTRACT 

Using the high-resolution spectrometer SPI on board the International Gamma-Ray Astro- 
physics Laboratory (INTEGRAL), we search for a spectral line produced by a dark matter 
(DM) particle with a mass in the range AQkcV Mjjm ^ 14A/eV^, decaying in the DM halo 
of the Milky Way. To distinguish the DM decay line from numerous instrumental lines found 
in the SPI background spectrum, we study the dependence of the intensity of the line signal 
on the offset of the SPI pointing from the direction toward the Galactic Centre. After a critical 
analysis of the uncertainties of the DM density profile in the inner Galaxy, we find that the 
intensity of the DM decay line should decrease by at least a factor of 3 when the offset from 
the Galactic Centre increases from 0° to 180°. We find that such a pronounced variation of 
the line flux across the sky is not observed for any line, detected with a significance higher 
than 3(7 in the SPI background spectrum. Possible DM decay origin is not ruled out only for 
the unidentified spectral lines, having low (~ 3a) significance or coinciding in position with 
the instrumental ones. In the energy interval from 20 keV to 7 MeV, we derive restrictions on 
the DM decay line flux, implied by the (non-)detection of the DM decay line. For a particular 
DM candidate, the sterile neutrino of mass Mdm, we derive a bound on the mixing angle. 
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1 INTRODUCTION 

Dark matter in the Universe 

There is a vast body of evidence, suggesting that the large frac- 
tion of matter in the Universe exists in the form of the Dark mat- 
ter (DM). However, while the total density of the D M is measured 
with a very high precision (Q.-oMii'^ = 0.105"*"° Qgl. lSpergel et al.l 



|2007|) . little is known about its properties apart from this. The pos- 
sibility that the DM is composed of the Standard Model (SM) par- 
ticles has been ruled out for a long time already. Indeed, the DM 
cannot be made out of baryons, as producing such an amount of 
baryonic matter would require drastic modifications of the sce- 
nario of the Big Bang nucleosynthesis (BBN), which otherwise 
successful ly describ es the abundance of light elements (see for 
example lDailll99^ . Recent microlensing experiments rule out 
the possibility that another type of baryonic DM - 
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pact halo objects (MACHOs) - con s titute dominant frac t ion of 
mass in the halo ( I Alcocketal] 120001 : iLasserre et al]|2000l ; lAlarj 
Il999l) . The only non-baryonic DM candidate in the SM can- 
didates - (left-handed) neutrino - is ruled o ut from the large 
scale structure (LSS) considerations (see e .g. iBond et aL 1 1 19801: 
iHannestad & Raffeltll2004 ICrottv et al.ll2004) . 

What are the properties of a successful DM candidate? 
First of all, this particle should be massive. Many extensions 
of the SM present the DM candidates with the masses rang- 
ing from ~ 10~^° eV (massive gravitons. iDubovskv et al. I l200^ 
and ~ 10"*^ eV (axions) to h undreds of GeV (WIMPs) and 
even to 10^=^ GeV (WIMPZILLA. lKuzmin & Tkachev l' 1998'. 19991: 
IChung etal.l 19991) For a r eview of particle phy s ics DM candidate s 
see e.g. lBergstroml ( l200dh : lBertone et"ai] ( l2005h : ICarr et alj ( l200^ . 

Secondly, there should exist mechanisms of DM production 
with the correct abundances. The production mechanism in partic- 
ular determines the velocity distribution of particles in the early 
Universe. This velocity distribution can, in principle, be probed ex- 
perimentally. Namely, if during the structure formation epoch the 
DM particles have velocities, comparable to the speed of sound 
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in the baryon-photon plasma, they "erase" density fluctuations at 
scales, smaller than the distance, they have traveled (called the 
free-streaming length). To differentiate various models in accor- 
dance with this property, the DM candidates with the negligible ve- 
locity dispersion (and, correspondingly, free-streaming) are called 
cold DM (CDM), while those with the free-streaming of the or- 
der of ~ 1 Mpc are considered to be warm (WDM)EI It is pos- 
sible to constrain the free-streaming length of a particular DM 
candidate by probing the structure of the Universe at galaxy-size 
scales. T his can be done through the analysis of the Lyman-a for- 



est data ( iHui et al.|[T997l) . Lyman-a analysis puts an upper bound 
on the free- stre a ming of the DM particles |Hansen et al] |2002|; 
Viel et alj l2005l : ISeliak et alj |2006| ; IViel et alj l2006l : IViel et atl 



2007h . It should be noted however that currently existing interpreta 



tion of the Lyman-a data is model-de pendent, as, apar t from a num- 
ber of astrophysical assumptions (see lHui et alj|l997h and compli- 
cated hydrodynamic simulations, it relies on a priori assumptions 
about the velocity distribution of the DM particles. 

A way to differentiate between CDM and WDM models 
would be to compare the numerical simulations of the DM dis- 
tribution in the Milky Way-type galaxies with the actual obser- 
vations. However, the resolution of the N-body simulations is not 
yet sufficient to answer the questions about e.g. the DM density 
profiles in dwarf satellite galaxies. Moreover, most of the simu- 
lations include only collisionless DM particles, and do not model 
the baryons and their feedback on the galaxy structure formation. 
These problems are not solved even for the CDM simulations, and 
WDM simulations have additional serious difficulties. From an ob- 
servational point of view, it has been argued for some time already 
that there is a discrepancy between CDM simu l ations and observa- 
tions (see e.g. Moore 1994:'Moore et al."l999:'Klvpin et al.''l999l; 
iBode et al. 2001 ; Avila-Reese et al.,2001 ; ,Goerdt et al..2006,) It has 
been claimed recently that a number of recent observations of dwarf 
satellite galaxies of the Milky way and Andromeda galaxy seem 
to indi cate the existence of the smallest scale at which the DM 
exists ( lGilmoreetid]|2006l . l2007l : iGilmord I2OO7I : IKodosov et al] 
l2007h . However, this statement and th e interpretation of the ob- 



servations are still subject to debate JKlimentowski et al 



IPenarrubia etal]|2007l ; Istrigari et al]|2007l : ISimon & GehjBoOTf) 



20071; 



Therefore it is too early to say what kind of DM models is favoured 
by comparing simulations and observations. 

Usually it is also necessary for the DM candidate to be stable. 
For the most popular DM candidate - weakly interacting massive 
particles (WIMPs), this is related to the fact that the particles of ~ 
electroweak mass, having weak strength interaction with SM matter 
(required to produce the correct amount of DM), would decay too 
fast and would not be "dark". If, however, the DM particle interacts 
with the SM more weakly than WIMPs, it could well have a finite 
(although cosmologically long) life time. 

There exis t several uns t able (decaying) DM can didates 
e.g. gravitino 1 Borgani et al. 19961; Isa ltz & Muravam^ I2OO3I; 



Roszkowski et al.l 



2OO5I; ICerdeno et airi2006i; ,Cembranos et al.l 



20061 ; iLola et alj|2007h . In this paper we will concentrate mainly 
on one candidate, the sterile neutrino (although our results 
will be applicable for any type of decay ing DM). Constraints 
on th e decaying DM w er e ana ly zed in [de Ruiula & Glashowl 



1 980l) ; iBerezhiani et id] ( Il987l); iDoroshkevich et al.l 
iBerezhiani et alj i 19901) ; Berezhiani & Khlopol 



198S); 



1990h 



iBertone et al. (l2007h ; IZhang et alj | |2007|) (see also the book 
bv lKhloDovlll997t) . 



Sterile neutrino DM 

It was noticed long ago that the right-handed (or as it is often 
called sterile) neutrino with the m ass in the keV range would rep- 
resent a viable DM candidate (Dodelson & Wid row 1994). Such 
a neutrino would interact with the rest of the matter only via the 
quadratic mixing with left-handed (active) neutrinos and therefore 
(although not stable) could have cosmologically long life-time. At 
the same time, it co uld be produc ed in the early Universe with the 
correct abundances foodelson & W idro w 1994; Shi & Fuller 199i; 
IShaposhnikov & Tkachevll 20061) . One of the decay channels of the 
unstable sterile neutrinos includes emission of photons of the en- 
ergy equal to half of the sterile neutrino rest energy. This potentially 
provides a possibility to observe the decays of DM sterile neutrinos 
via detection of a characteristic spectral line in the spectra of astro- 
physical objects with large DM concentration. 

Recently this DM ca ndidate has attracted much attention (see 
e.g. lShaposhnikovl klWlt ) and references therein). It was found that 
a very modest and natural extension of the SM by 3 right-handed 
neutrinos (making the SM more symmetric as all SM fermions, in- 
cluding neutrino, would have now their left and right handed coun- 
terparts) provided a viable extension of the theory, capable of solv- 
ing several "beyond the SM" problems. First of all, such an exten- 
sion makes neutrinos massive and thus perhaps provides the sim- 
plest and the most natural e xplanation of the p h enomenon of "neu- 
trino o s cillations" (see e.g. iFogli et alj ( l2006h ; IStrumia & Vissanil 
( I2OO6I) ; iGiuntil l l2007h for reviews). The smallness of neutrino 
masses in this model (called z/MSM in Asaka & Shaposhniko^ 
I2OO5I) is achieved by the usual see-saw mechanism with Majorana 
masses of right-handed neutrinos being below electroweak scale0 

Secondly, if two heavier sterile neutrinos {N2 and N^) 
are almost degenerate in mass and have their masses between 
O(IOO) MeV and 0(20) GeV, the uMSM provides the mechanism 
of generating the baryon asymmetry of the Universe. Thirdly, the 
lightest sterile neutrino TVi can have arbitrary mass and arbitrarily 
weak coupling with the (active) neutrino sector. At the same time, 
it can be produced in the early Universe in the correct amounts. It 
represents therefore the DM particle in the z^MSM. Thus, altogether 
the vMSM represents (arguably) the simplest extension of the SM, 
capable of explaining three important questions: origin and small- 
ness of neutrino masses, baryon asymmetry in the Universe and the 
existence of the DM. 



Existing restrictions on sterile neutrino DM parameters. 

What are the current restrictions on parameters (mass and mixing) 
of sterile neutrino DM? First of all sterile neutrino mass should 



^ The left-handed neutrino would represent hot DM in this terminology, 
i.e. the DM with the free-streaming length ;§> 1 Mpc. 



^ The fact that the i/MSM does not introduce any new scale above the 
electroweak one, makes this theory especially appealing from the point of 
view of its experimental verification/falsification. 
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satisfy the universal Tremaine-Gunn lower bound[f| Mdm }t 300 — 
500 eVQ 

Next, as the sterile neutrino possesses the (two-body) radiative 
decay channel: A^i — > ly+j, the emitted photon would carry the en- 
ergy E-y = Mdm/2. a large flux of such photons is expected from 
the large concentrations of the DM sterile neutrinos, like galaxies 
or galaxy clusters. 

Recently an extensive search of the DM decay line in the 
region of masses Mdm 5, 20 keV was conducted, using the 



data of Chandra llRiemer-S0rensen et a D l2006l : 
2006d|: lAbazaiian et al.ll2007h and XMM-Newton 
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2006allbllct IWatson et alj |2006| ; ISovarskv et al] l2007h . The re- 
gi on of soft X-ray (dow n to energies 0.2 keV) was explored 
by iBovarskv et alj 1 2007 ) with the use of the wide field of view 
spectrometer ( McCammon et al.l l2002h . The non-observation of 
the DM decay line in X-ray, combined with the fir st principles 
calculation of DM prod uction in the early Un i verse jAsaka et al.l 
l2007h . implies that the iDodelson & Widrowl ( 1 19941) (DW) sce- 
nario c an work only if th e sterile neutrino mass is below 
4 keV teovarskv et al.l l2008h . If one takes into account recent 
lower bound on the mas s of s terile neutrino DM in the DW sce- 
nario A/dm ^ 5.6 keV jViel et al. 2007), it seems that the pos- 
sibil ity that all the DM is produced via DW scenario is ruled 
out teovarskv et alj|2008l) . The possibility that only fraction of the 
DM i s produced via DW mechanism remains open jPalazzo et al.l 
l2007h . 

There are other viable mechanisms of DM production, in- 
cluding e.g. reson ant oscillation prod uction in the presence of lep- 
ton asymmetries ( IShi&Fulleilll999h . Sterile neutrino DM can be 
produ ced by the decay of light inflaton jShaposhnikov & TkachevI 
l2006l) or in a similar model with th e different choice of parame- 
ters ( lKusenkdl2006l J Petraki & Kusenko 2007) . These mechanisms 
are currently not constrained and remain valid for DM particles 
with the masses in the keV range and above. 

The search for the DM decay line signal produced by ster- 
ile neutrinos with masses above ~ 20 keV is complicated by the 
absence of the focusing optics telescopes (similar to Chandra or 
XMM-Newton) in the hard X-ray and 7-ray domain of the spec- 
trum. For ex ample, the ex isting restrictions in the 20 — 100 keV 
mass range teovarskv et a l. 2006a are derived from the obser- 
vations of diffuse X-r ay background, with the help of non-imaging 
instruments, HEAO-I l lGniberetal.ll999l) . The current status of as- 
trophysical observations in summarized in lRuchavskivl ( l2007h . 

In this paper we use the spectrometer SPI on board of IN- 
TEGRAL satellite to place restrictions on parameters of decay- 
ing DM in the mass range 40 keV — 14 MeV. This range 
of masses is interesting, for example, the sterile neutrinos, pro- 
duced in the early Universe in the presence of large lepton 

J ^ 1 I 1 ^ I 

asym metries fShi & Fuller 1999) or through the inflaton de- 
cay dShaposhnikoy & Tkachev 200^. It is also relevarit for t he case 
of gravitino DM [Pagels & PrimacJl982l : lBond et aljl982h . 
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^ In its simplest form the Tremaine-Gunn bound comes from the 
fact t hat for the fermions th e re is a maxi mal density in the phase 
space iXremaine & Gunnlll979t rPalcanton & H ogan 2001) and therefore 
the observed phase-space density in various DM dominated systems should 

be less that this (mass dependent) bound; 

A stronger lo wer bound from Ly-a iSehak et al]|2006l : Iviel et alj|2006l : 
IViel et al.l I2OO7I) can be obtained in the case of t he particular produc- 
tion m echanisms - the Dodelson-Widrow scenario toodelson & Widrow' 
I994I). For other possible produ ction mechanisms (e.g. Shi & Fuller 1999; 
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Shaposhnikov & Tkachevl2006l) the Ly-Q constraints should be reanalyzed. 
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Figure 1. Comparison of sensitivity towards the search of the narrow DM 
decay line for different instiuments with the wide FoV. Diagonal straight 
lines show the improvement of sensitivity (by a factor, marked on the line) 
as compared with the HEAO-I A4 low energy detector (LED), taken as a 
reference. 



\ yhen the pr e parati on of this paper was at its final 
stage, lYuksel et al.l ( |2008 | , hereafter Y07) published th eir work, 
which used the results of iTeegarden & Watanabd ( I2OO6I hereafter 
T06) to place restrictions on the parameters of sterile neutrino DM 
in the range 40 — 700 keV. We discuss it in more details in Sec- 
tion|6] 



SPI spectrometer 

The absence of the focusing optics significantly reduces the sensi- 
tivity of the telescopes operating in the hard X-ray/soft 7-ray en- 
ergy band. Most of the instruments operating in this energy band 
use collimators and/or coded masks to distinguish signals from the 
sources on the sky from the instrumental background. Contrary to 
the focusing optics telescopes, both the source and background sig- 
nals are collected from the entire detector, which significantly in- 
creases the irreducible background. 

The focusing optics enables to significantly reduce the back- 
ground only in the studies of point sources. If the source under in- 
vestigation occupies a large fraction of the sky (e.g. the entire Milky 
Way galaxy), the performance of the focusing and non-focusing in- 
struments with the same detector collection area are, in fact, com- 
parable. 

In the case of an extended source, emitting a narrow spec- 
tral line, an efficient way of reduction of instrumental background 
is via the improvement of the spectral resolution of the instru- 
ment (in the case of a broad continuum background spectrum, the 
number of background counts at the energy of the line is propor- 
tional to the spectral resolution AE). The best possible sensitivity 
is achieved when the spectral resolution reaches the intrinsic width 
of t he spectral line ( see Fi g[T]for the case of wide FoV instruments 
and lBovarskv et aL ||2003) for the case of narrow FoV instruments). 

In the case of the line produced by the DM decaying in the 
Milky Way halo, the line width is determined by the Doppler broad- 
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Figure 2. The geometry of the SPI FoV. 



ening by the random motion of the DM particles. The velocity dis- 
persion of the DM motion in the halo is about the rotation velocity 
of the Galactic disk, v ~ 200 km/s. This means that Doppler broad- 
ening of the DM decay line is about 



V Q 



(1) 



Thus, the optimal spectral resolution of an instrument searching for 
the DM decay line produced by the Milky Way DM halo should be 

Such optimal spectral resolution is almost achieved with the 
spectrometer SPI on board of INTEGRAL satellite, which has the 
maximal spectral resolving power of E/AE c::; 500 and w orks 
in the energy range 20 keV - 8 MeV jVedrenne etal]|2003l) . SPI 
is a "coded mask" type instrument with an array of 19 hexagonal 
shaped Ge detectors (of which only 17 are operating at the mo- 
ment). 

The SPI telescope consists of a coded mask inscribed into a 
circle of the radius i?mask = 39 cm, placed at the height H = 
171 cm above the detector plane and of the detector, which has the 
shape of a hexagon inscribed into a circle of the radius i?dct — 
15.3 cm (see Fig.|2j. The portion of the sky visible from each point 
of the SPI detector (the so-called /ifZ/y coded field of view, FCFOV) 
has therefore angular diameter 



©FCFi 



2 arctan 



R„ 



Rd 



H 



16° 



(2) 



while the portion of the sky visible by at least some of the detectors 
(the partially coded field of view, PCFOV) is 



Op 



= 2 arctan 



Rn 



+ Rd 



H 



:35° 



(3) 



The solid angle spanned by the cone with this opening angle is 
rJpcFov = 2tv(i- cos(epcFov/2) j ~ 0.29 (see Fig.|2ll- Wide field 
of view makes the SPI telescope suitable for the study of the very 
extended sources, like the Milky Way DM halo. 




Figure 3. The effective area of the SPI detector for an on-axis source, as a 
function of the photon energy. The plot is produced by collective the on-axis 
effective areas of the 17 SPI detectors from the instrumental characteristics 
files. 



2 THE EXPECTED SIGNAL FROM THE DM DECAY IN 
THE HALO OF THE MILKY WAY. 

The expected surface brightness of the DM decay line in a given 
direction on the sky is a function of the angular distance between 
the given direction on the sky and the direction towards the Galactic 
center (GC). It can be calculated by taking the integral of the DM 
density profile Pdm('") along the line of sight ("column density") 



5dm {(f>) = dzp, 



■ 2zrQ cos <j) -\- z'' 



(4) 



where tq ^ 
GC . Angle 



8.5 kpc is the distance from the Solar system to the 
b is related to the galactic coordinates (b, I) via 

cos S = cos b cos I . 



(5) 

Thus, the galactic center corresponds to = 0°, the anti-center 
(f> — 180°, and the direction perpendicular to the galactic plane to 
— 90°. The expected DM flux is given then by 

dFoM (0) _ FdmSt, 



-Sdm{<P) 



(6) 



dn 47rAfDM 
where Fdm is the DM decay rate 

In general, the surface brightness Fdm(0) is variable across 
the telescope FoV. This is especially true for a wide field of view 
(FoV) instruments (like SPI). In order to calculate the detector 
count rate, one has to integrate flux ^ over the FoV and over the 
(effective) detector area and then divide by the energy of the pho- 
tons, Ej = Mdm/2: 



dad [3 



At,a{Ej\a, 13) dFau 



E~^ 



dVt 



{4>{a,P)) , (7) 



where (q, /3) are the angular coordinates in the FoV, Acff is the 
effective area at energy E^ for the photons, coming from the direc- 
tion (a, /3). 

The effective area of the SPI detector (which is determined by 
the transparency of the mask and the quantum efficiency of the de- 
tector) changes with the photon energy. For an on-axis point source, 

dF/dn{a, 13) = foS{a)S{f3) , 

the integral of Eq. ^ reduces to /oj4cff,on, where Acd .on{E.y) is 
the detector effective area/or an on-axis source. Its dependence on 
energy Ej is shown on Fig. [30 

^ The on-axis effective area is calculated by summing the energy- 
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Figure 4. Dependence of the effective area on the off-axis position of a 
(point) source. 



In the general case of extended sources, evaluation of tlie de- 
tector count rate (O analytically is not possible because of the 
complicated dependence of the effective area on the off-axis angle 
(shown on Fig.|4ll. In the simplest case of an extended source with 
a constant surface brightness dFomift') /d^l = /ext = const, the 
integral of Eq. Q reduces to the multiplication by the solid angle 
ilpcFov — 0.29 and the effective area, averaged over the FoV: 

^eff,ext(-B^) = -TT^ — / [ dad/S Aaf{{E^\a, f3) 

' ipCFOV J J 



FoV 

K{Ej)Acff,ou(Ej) 



(8) 



The numerical factor k{E^) depends on the energy and has to be 
calculated via a numerical integration over the energy dependent 
off-axis response map of the SPI detector. A reasonably accurate 
numerical approximation to k,[E^) is given by 



fii0.165(£/keV)" 



(9) 



One can see that k ^ 1 in all the energy interval. This is explained 
by the fact that the detector area visible from a given direction on 
the sky strongly decreases with the increase of the off-axis angle 
of this direction, so that the sky-averaged effective area is much 
smaller than the on-axis effective area of the detector. Substituting 
([9]l, {Sjl into ([7} one finds that for an extended source of constant 
surface brightness the detector count rate is 



Re: 



2.73 X 10"' 

^cff,exl(-B7 



cts 

S 



IkeV 



150 cm2 



(dFr,M/dn), 



(10) 



10~^^erg/(cm2 ssr) 



dependent on-axis effective areas of each of the 17 operating detectors of 
SPI, extracted from the instrument's characteristics files. 



2.1 Modeling the DM halo of the Galaxy 

The DM halo of the Galaxy has been extensively studied (se e e.g. 
iKravtsov et al .ll998l ; lKlvDin et alj2002l ; lBattaglia et alj2005l) . Var- 
ious DM profiles, used to fit observed velocity distributions, differ 

the most in the GC re gion. 

It was shown in iKlvpin et~al] ( l2002h : iBattagUa et"al] ( |2005|) 
that the DM halo of the M W can be described by the Navarro- 
Frenk-White (NFW) profile jNavarro et al]|l997h 



Pnfw (?") 



r(r + Ts 



(11) 



with parameters, given in Table [T]The relation between virial pa- 
ra meters and ps, Va can easily be found (see e.g. the Appendix A 
of lBovarskv et al.ll2007h . 

To explore the uncertainty of the DM density profile in the 
inner part of the Galaxy, we also describe the DM distribu tion in 
the MW via an isothermal profile ( lBahcall&Soneiralll980l) : 



Po 



AuGn + r1 



1 + {rlr^Y 



(12) 



The following parameters of isothermal profile reproduce the DM 
contribution to the (outer part s of) G alaxy rotation curve = 
170 km/sec and = 4 kpc l lBovarskv e t al. 2006c, 200^) (i.e. 

po = 1.2 X 10' 



,6 keV 
cm' 



_j j^ 4kpc j ^ xjjese parameters are 



1^170 km/s J L - J 

consis t ent with those, f rom f avored NFW models of iKlvpin et al.l 
ll2002h : lBattaglia et alj (l2005h . i.e. for 4> ^ 90° the difference be- 
tween isothermal model and NFW with preferred parameters was 
completely negligible (less than 5%) - c.f. F1G[5] Both types of 
models provide the local DM density at the position of the Sun to 
be Pdm (7'f)) — 0.22 GeV/cm'^, w hic h is close to the existing es ti- 
mates l lKuiiken & Gilmorelll989cllibl Il99ll : iGilmore et aljl989l) . 

The DM flux from a given direction 0, measured by an ob- 
server on Earth (distance tq ~ 8.5 kpc from the GC), is given 
by 



>Siso(0) ■ 



R 



+ arctan 



arctan 



/ rQ cos tt> \ 

( « ^ 

rQ I cos 4>\ J ' 



COS (j) 
COS (f) < 



(13) 

where R = \Jr'l + r'^ sin^ (j> and p^rc 1.5 x 10^*^ keV/ cm^. 

The uncertainty of the DM radial density profile in the inner 
Galaxy stems from the difficulty of separation between visible and 
DM contributions to the inner Galaxy rotation curved In order to 
get the most conservative limit on the column density of the DM in 
the direction of the GC, one can assume the following "rigid lower 
bound"; while the DM outsi de the rp, is describ ed by the "maximal 
disk" model (model A2 of IKlvpin et"aDl2002h . for r ^ r© DM 
density remains constant (so t hat the total DM m ass within r© is 
the same as in the model A2 o f lKlvpin et al]|2002l) . This gives 

^6 Mq 



mm 



3.9 x 10° 



kpc^' 



: 0.146 X 10' 



keV 



(14) 



The surface brightness profile on the "constant density" model 
is shown in black dashed line on the Fig|5] One can see that 
the difference between the maximal ((f) = 0°) and the minimal 
(j) = 180°) column densities is ~ 3.4 (as compared to ~ 6 for 



^ When quoting results of lKlvpin et alj )2002h . we do not take the effects 
of baryon compression on DM into account. While these effects make DM 
distribution in the core of the MW denser, any such computation is strongly 
model dependent. 
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Table 1. Best-fit parameters of NFW model of the MW DM halo. Max. disk model maximizes amount of baryonic matter in the inner 3 kpc of the MW halo 

(A/DM/(Mdi,sk + Afbulgc) = 0.4 for the model A2 and J\/DM/(Mdi.sk + A4i,lgc) = 0.14 in the model B2). 



References 


A/vir [A^0] 


rvir [kpc] 


Concentration 


Ts [kpc] 


Ps [A/0/kpc^] 


Klvpin et al. (2002). favored models (Ax or Bx) 


1.0 X 10^^ 


258 


12 


21.5 


4.9 X 10" 


KlvDin et al. (2002), Max. disk models Ao 


0.71 X 10^2 


230 


5 


46 


0.6 X 10'' 


Klvpin et al. (2002), Max. disk models Bo 


0.71 X 10^2 


230 


10 


23 


3.1 X lO'' 


BattagUa et al. (2005) 


O.Slg 2 X 10^2 


255 


18 


14.2 


11.2 X lO'' 




100 
^[degreesj 

Figure 5. Expected column density for vaiious DM profiles: favored NFW 
profile (red thick solid line); NFW profile with the maximal disk (model 
A2, see Table[T) - blue solid line; cored (isothermal) profile - green thick 
dashed line; constant density within tq - black dashed line). 



isothermal model). For comparison we show on Fig. [5] expected 
DM flux © for various profiles. The minimal column density is 
of course the one in the direction of anti-center: S[<j) = 180°) ~ 
0.33 X 10^® keV/ cm . We see that even for the minimal profile 

S((/.<30°) ^ lO^^^keV/ cm^ 



2.2 DIM decay line count rate 

In the case of the Majorana sterile neutrinos of mass AfpM the 
DM d ecay width is given by l iPal & Wolfensteial982 : iBarger et al.l 



Fdm ==: 1.3 X 10" 

Substituting dlSb to ^ we find 

dFotA 



sin^ 261 



10- 



A^DI 



IkeV 



(15) 



dQ. 



3 X 10 erg/(cm^ s sr) 



(16) 



sin^ 261 



10- 



IkeV 



5dm (0) 



1028 keV/cm^ 



The quoted value of Fdm is for the Majorana sterile neutrino. In ca se of 
Dirac particle this value is 2 times smaller (c.f. |Pal & Wolfenstein 19821 : 
lBargeretal.ll995l) . 



The lower bound on the DM decay line rate in SPI pointings toward 
the inner Galaxy is calculated by substituting the column density 
5 = 10^** keV/ cm^ (see Fig.|5j into Eqs. ([lOj, ([16) 



3.0 X 10" 



cts 



5dm (<?!') 



10^8 keV/cm^ 



sin^ 261 



10- 



1 keV 
(17) 

The approximation of the constant surface brightness works well, 
if the extended source has a core of the angular diameter exceeding 
the size of the SPI partially coded FoV (Spcfov ~ 17° maximal 
off-axis angle). Taking isothermal profile the angular size of the flat 
core of the extended source is 



(pcore = arctan(r'c/rQ) ~ 25° , 
which satisfies this constraint. 



(18) 



3 STRATEGY OF SEARCH FOR THE DM DECAY LINE 
WITH SPI 

The MW halo contribution to the DM decay signal represents the 
all-sky source. Indeed, as the results of Section \2A] show, the vari- 
ability of the signal over the sky may be as low as the factor ~ 3. 
This makes the strategy of search of the DM decay signal different 
from any other types of astrophysical sources: the point sources, 
diffuse sources (e.g. ~ 10° Gau ssian profile for e'^e~ annihilation 
region, iKnodlseder et alj dlOOSh) or even the search for DM anni- 



hilation signal (see e .g.'Tasitsiomi et al. 2004; Boehm et al. '2004|; 



iDiemand 



Lgnal (see e .g. 
et al.l2007l :ISi 



Sanchez-Conde e t al. 200_6; Carr et al. 2006|) . 



The problem gets exacerbated by the fact that during its mo- 
tion, SPI is irradiated by the charged high-energy particles (parti- 
cles from Earth radiation belt. Solar wind, cosmic TeV photons). As 
a result, the materials (even detectors themselves) used for SPI con- 
struction start to radiate in different energy regions (see subsection 
13.21 ). As a result any SPI spectrum consists of a broad continuum, 
which is a combination of the sky and instrumental backgrounds, 
and of a set of the instrumental background lines ( Attie et al. 200^ 
iDiehl et alj2"003l : |jean et al.l2003l : IWeidenspointner et al.l2003l) . In 
order to detect a spectral line produced by an astrophysical source 
one has to be able to (a) separate the continuum and line contri- 
butions to the spectrum and (b) separate the instrumental and sky 
signal contributions to the lines found. 

One can expect three a priori situations: 

(I) DM decay line is strong (its equivalent width much larger 
than the spectral resolution) and at its position there are no 
other strong lines (of either instrumental or astrophysical ori- 
gin). Such a line, due to its presence in any SPI spectrum and 
its low variability over the sky can in principle be confused 
with some unknown instrumental line. 
(II) DM line is weak (~ 3 — 4(t detection over the continuum) 
but its position also does not coincide with any instrumental 
line. 



Constraining DM properties with SPI 7 



(III) DM decay line coincides with some instrumental line. To 
be able find such a line we need to model SPI instrumental 
background. 

To be able to work effectively with all these situations, we need to 
find the way to separate the source and background contributions. 



3.1 Imaging 

To distinguish source and background contribution to the signal, 
one often uses imaging capabilities of an instrument. If the size 
of a point or even an extended source on the sky is smaller than 
the size of the SPI FoV one can (at least, to some extent) use 
the imaging capabilities of the SPI instrument. In this case the 
coded mask, placed above the detector, partially screens the in- 
dividual detectors from the source, so that the source at a given 
position on the sky produces different count rates in different de- 
tectors. One can find the source flux by comparing the ratios of 
the actual count rates in different modules of the detector to the 
ones predi cted by the degree of screening of the module s by the 
mask (see iDubath et all |2005| ; ISkinner & Connelll l2003h . It is a 
challenge, however, to use the imaging capabilities of the SPI to 
separate the astrophysical signal from the instrumental background 
if the size of the extended source is comparab le to the size of the 
SPI FoV (see e.g. iKnodlseder et al.l l2005: A Uain & Roauesll2006l ; 
IWeidenspointner et al.l 2007 . and refs. therein). Therefore for our 
analysis we did not use any imaging capabilities of SPI, and to pro- 
duce spectra from some point in the sky we just collected all the 
photons, arriving in the SPI FoV. 



3.2 SPI background modeling 

In the absence of imaging, the separation of the instrumental and 
astrophysical contributions to the line spectrum req uires some sort 
of background mode l ing (see e.g. Weidenspointn er et al.l 1200 3l ; 
iTeegarden et al.ll2004l : iTeegarden &'WatanabQi2006i) . Namely, for 
the background modeling we can use the fact that for any DM dis- 
tribution model the intensity of the DM decay line changes by a fac- 
tor ^ 3 between the pointings towards the Galactic center (0 ~ 0°) 
and anti-center ((j> ~ 180°, see Sec |2.1| (. On the other hand, if the 
line is of purely instrumental origin, there is no a-priori reason why 
the strength of the line in the background spectra of the pointings 
towards e.g. the Galactic Anti-center should be different from the 
strength of the line in the spectra of the pointings toward e.g the 
GC. Thus, one possible way to distinguish between the DM decay 
and instrumental origin of the line is to study the variations of the 
line's strength depending on its sky position (in the simplest case - 
on the "off-GC" angle (j>, of the pointing, Eq. l|5j)- 

The situation becomes more complicated due to the fact that 
the instrumental background (and thus the intensity of the in- 
strumental lines) experiences great variability in time (depend- 
ing on the position in orbit, solar flares and the solar activity 
period, degradation o f the detectors, etc., c.f. Ijean et al] 1200 3l ; 
ITeegarden et al.|[2004h . As observations of different parts of the 
sky can be significantly separated in time, one needs to use "back- 
ground tracer s" to find the correct spatial dependence of the line 
intensity dJean et a l. 2003; Teegarden & Watanabe 2006). Without 
some sort of "renormalization" procedure, which corrects the ab- 
solute value of the line flux using a measurement of a specific 
characteristics of the SPI instrument as a "calibrator" of the flux, 
the (j!> dependence for any of the detected lines contains no useful 



information. There exist various "background tracers" (Ge detec- 
tors saturation rates, a nti-coincidence sh i eld rates, rates of certai n 
background lines , see 'Jean et alJ j2003l) : ITeegarden et"ai] ( |2004|) : 
ITeegarden & Watanabe. ^2006) and refs. therein). 



3.3 Searching for the lines 

To be able to detect strong DM line, which is not close in position 
to any instrumental line (case I above), we used the modification of 
the method of background subtraction, described in TW06. TW06 
looked for 7-ray lines, assuming different types of sources, from 
the point sources to the very diffuse sources (10° Gaussian, 30° 
flat, etc.) TW06 showed that the strong background line at 198 keV 
can be used as a background tracer, if background observations are 
matched close in time to the corresponding "source" ones. This al- 
lowed TW06 to cancel all strong instrumental lines with the preci- 
sion better than 1%. TW06 detected no emission line in such back- 
ground subtracted spectrum (apart from the 511 keV and 1 809 keV) 
with the significance above 3.5a. 

We adopt the following modification of the TW06 method: 
- As the DM decay signal remains nearly constant within cen- 
tral 30 — 50° , the method of TW06, if applied directly, could cancel 
most of the DM signallf] We therefore subtract the data (renormal- 
ized by the strength of 198 keV line) in the direction away from the 
GC (off-GC angle > 120°) from the ON-GC dataset (the angle 
13°). 

~ In the resulting "ON-OFF" spectrum we perform the search 
for the line with the significance higher than 3cr. 

This procedure allows to eliminate strong instrumental lines 
with the precision better than few percents. At the same time any 
strong DM line would remain in the "ON-OFF" spectrum. Indeed, 
even for the flattest profile (Section |2. lb , the strength of the DM 
signal in the OFF dataset is at least 60% weaker that of the ON 
dataset. Therefore we see, that the modification, described above, 
is indeed well suited for searching of the strong DM decay line 
(case I). 

However, this method does not work well for the weak (3— 4cr) 
lines, or for the lines, whose position coincides with some instru- 
mental line (cases II-III above). Indeed, in this case it is not pos- 
sible to tell whether the remaining line is the residual of the in- 
strumental one or has the astrophysical origin. Below we will use 
an alternative method of analysis of the detected lines, suitable for 
cases II and III. 



3.4 Analyzing a candidate line 

Having detected a number of lines with the significance of 3(7 and 
above, we should decide which of them can be considered as "DM 
decay line candidates". To this end we do the following. 

a) We compare line flux for each of these lines with the flux 
of the same line in the "ON" spectrum. We decide that the line is 
a "DM line candidate" if the cancellation of the flux between ON 
and OFF datasets was worse than 1O%0 



° For example for the most conservative DM distribution model, the differ- 
ence of DM signals at = 0° and </> = 30° is mere 8%. 
^ In principle, the DM line in ON-OFF spectrum should not cancel by 
more than ~ 40%, while the background instrumental line should cancel 
better than 1%. Thus the choice of the threshold to be ai'ound 10% ensures 
that no DM decay line was thrown away while most of the instrumental 
lines disappeared. 
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b) For any "DM candidate line" we construct its "spatial pro- 
file" (as described in details in the next Section) to check for the 
possibility of it to be a DM decay line (we also construct distri- 
butio n of the line flux over the sky for all the unidentified lines 
from lWeidenspointner et al .l2003h . Since the column density of the 
DM in the direction toward the GC should be higher than that of in 
the direction toward the Galactic anti-center, one should see a grad- 
ual decrease of the line strength with the increasing angle 4>- We do 
not make any specific assumption about the DM density profile and 
do not try to fit the candidate line spatial profile to any particular 
model, but rather look if there is a general trend of decreasing in- 
tensity of the line with the increasing off-GC angle. 



4 DATA REDUCTION 
4.1 ON dataset 

During its almost 5 years in orbit INTEGRAL has intensively ob- 
served the inner part of the Galaxy (Galactic Center, Galactic 
Bulge and the iimer part of the Galactic Plane) and collected about 
Texp ~ 10 Ms of exposure time in the GC region. In our analysis 
of the inner Galaxy we used the publicly available data (as of July 
2007) from all INTEGRAL pointings at which the angle off the GC 
was at most 13° and for which the SPI exposure time was larger 
than 1 ksec. This criteria selects 5355 pointings (or "Science Win- 
dows", ScW), with total exposure time of 12.2 Ms, spread over the 
period from February 2003, till April, 25, 2006. We call this dataset 
"ON" dataset. 

For each of the analyzed ScWs, we have extracted photon 
(event) lists from spi-oper . f its files and applied additional 
energy correction to convert the channel number into photon en- 
ergy, using spi_gain_cor tool from standard Offline Analysis 
Software (OS A). We have binned the events into narrow energy 
bins of the size ASbin = 0.5 keV to generate the background 
counts spectra in each ScW, each revolution and, subsequently, in 
the entire data set. 

We then applied the "sliding spectral window" method to the 
line search (as described e.g. in TW06 to produce a continuum sub- 
tracted spectrum of the " ON " dataset. Namely, at each given energy 
Eq, one defines an energy interval Eo — 2AE < E < Eo + 2AE, 
where AE is the SPI spectral resolution at a given energy, as a 
"line signal" energy band. For the (energy dependent) AE we used 
approxi mate formula from SPI/INTEGRAL ground calibration of 
FWHM jAttie et al.l2003h : 



AE{E) = Fi+ F'zVe + F3E 



(19) 



where Fi = 1.54, F2 = 4.6 • 10"^, F3 = 6.0 • 10"" and energy 
E is in keV. For E = 10^ keV FWHM « 2.3 keV. 

For each energy bin centered at an energy Eq we have defined 
the two adjacent energy intervals, Eq — 4A_E < E < Eq — 2AE 
and Eo + 2A_E < E < Eo + 4A_E, and postulated that the sum of 
the count rates in these two adjacent energy bands gives the mea- 
sure of the continuum count rate in the energy band around Eq. 
Subtracting the sum of the count rates in the adjacent energy bands 
from the count rate in the "line signal" energy band, we have cal- 
culated the continuum subtracted count rate at a given energy Eg. 
Doing such procedure at all the energies 20 keV < -Bo < 8 MeV, 
we have produced a "continuum subtracted" SPI background spec- 
trum.In this spectr um we were able to identify most of the known 
instrumental lines jWeidenspointner et al .l2003h . 
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Figure 6. Comparison of ON-OFF spectrum (thick solid line) with the 
0.01 X the ON spectrum after the sliding window (thin dashed line). It 
can be seen that the instrumental lines are subtracted with the precision bet- 
ter than 1%. 



4.2 ON-OFF dataset 

Most of the lines found in the continuum subtracted background 
spectrum are of the instrumental origin. To remove them, we 
matched each ScW in the ON dataset with the pointing away from 
the GC (galactic coordinate (ji > 120°) - OFF pointing. As de- 
scribed by TW06, the 198 keV line can serve as good background 
tracer if the time duration between ON and OFF observations is 
< 20 days. We were able to match 3688 ON-OFF pairs. For 
each ON-OFF pair we introduced normalizing coefficient n for the 
OFF spectrum in such a way that the strong instrumental line at 
198 keV cancels completely after subtraction of the OFF spectrum 
multiplied by the factor n from the ON spectrum. After that we 
subtracted (renormalized) OFF ScW from the corresponding ScW 
from the ON dataset. This allowed us to remove most prominent 
instrumental lines with the precision better than 1% (c.f. Fig. |6]l. 
To avoid contributions of strong astrophysical sources (such as 
e.g. Crab) we threw out all pairs with negative total flux at 20- 
40 keV range after subtraction. Taking average over 2456 remain- 
ing "good" pairs we received the spectrum almost free from back- 
ground at energies above 200 keV. At low (< 200keV) energies we 
found continuum component, which can be fitted with the simple 
power law: 



FiE) = Fc 



E 



100 keV 



Parameters of this background were found to be 

Fo = (4.95 ± 0.05) X lO^-'^cts/s/cm^ 
a = -(2.264 ±0.003) 



(20) 



(21) 



This continuum represents the residual contribution from all the set 
of the astrophysical sources present in the Galactic Bulge. 

4.3 Systematic error 

To estimate the systematic error of our "ON-minus-OFF" dataset, 
we computed background around the "tracer line" of 198 keV. We 
found that it does not vanish. Thus, we estimated the systematic 
error as the error in the normalization coefficient n which would 
make the background zero within systematic uncertainty. This cor- 
rection Sn can be found as follows. Let n be the coefficient, needed 



Constraining DM properties with SPI 9 



Phase 


Revolutions (start-stop) 


1 


042-092 




096-140''^^^^^| 


2 


140-205 




209-215*^^^^^! 


3 


215-277 


4 


282-326 


5 


330-395 


6 


400-446 



Table 2. Splitting revolutions into phases in con'espondence with annealing 
phases and breakage of the detectors (revolutions, marked with the *). 



to cancel flux in the 198 keV line in ON and OFF spectra: 

n — Fon/Foff 



(22) 



where Fon,Foff- fluxes in 198keV line in ON and OFF ScWs 
correspondingly. The remaining non-zero SF flux in the adjacent 
to line position in ON-OFF spectrum, determines the uncertainty 
of the coefficient: 



5n ■ 



5F 

FoN 



(23) 



We found, that average value of (Sn) is equal to (Sn) = 1.1 • 
10"'', (n) ~ 1. So, our systematic error of final ON-OFF spectra 
atenergy Sis 1.1- 10^=^ Foff(-B) ^ 1.1 • IO^'^Foat. We add this 
systematic uncertainty to the flux of ON~OFF spectrum in every 
energy bin. 



4.4 Obtaining 3a restrictions 

At the energies at which no lines were detected (i.e. the "continuum 
subtracted" count rate did not deviate by more than 3a from zero) 
we obtained the 3a upper limit on the possible flux from the DM 
decay. Above ~ 200 keV the flux in the energy bin is zero within 
statistical errors, therefore 3a upper limit flux is given by statistical 
plus systematic errors. Below 200keV we put statistical restrictions 
above power law continuum flux i20\ . described in the Section l4~2l 
Using Eq. J17t one can derive the restriction on the sterile neutrino 
mixing angle, implied by this upper limit. One should also take into 
account that the subtraction of the OFF observations led to the re- 
duction of the expected DM signal. Taking the most conservative 
"minimal" model, described in Section lZT] we see that the subtrac- 
tion of the OFF signal leads to about 40% decrease of the expected 
DM signal0 The resulting 3a bound is shown on FigllOl 



4.5 Possible DM candidates 

When analyzing ON-OFF spectrum, we found that almost all lines, 
present in ON spectrum cancel with precision better than few per- 
cents. We found 21 lines (see Table 3) that did not cancel by at 
least 90%, (including known lines at 51 1 keV and 1809 keV). Apart 
from these 2 lines all other lines are detected with low significance 
3-4cr. 

As discussed in Sections l3.3H3.4l we took all these lines as 
possible DM candidates and analyzed the dependence of the line 
fluxes F{(j)) on the off-GC angle cj) of the pointing. If the DM dis- 
tribution in the inner part of the Galaxy were known, it would be 



' To estimate this, we took the maximal column density for OFF observa- 
tions ate* = 120° - 17°. 




600 



Figure 7. Position on the sky as a function of revolution over 6 years of 
INTEGRAL observations. The periods of annealing phases are shown in 
solid vertical lines. Two dashed lines indicate the revolution, duiing which 
2ND and 17th of 19 SPI detectors have failed. 
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Figure 8. Scatter of the flux data points for the line al E = 1068 kcV as a 
function of the off-GC angle. 



possible to distinguish between the instrumental and DM decay ori- 
gin of a line by fitting F{(j)) with a known profile calculated from 
the radial DM density profile. However, the details of the radial 
DM density profile in the inner Galaxy are highly uncertain, and 
this prevents us from directly fitting the model profile to the data. 
We adopted a simple criterion which selects a DM decay candidate 
line: the ratio of fluxes 



7^ = 



F(0°) 
F(180°; 



3. 



(24) 



where TZmm is the ratio of the DM decay line fluxes from the GC 
and the Galactic anti-center in the "minimal DM content" model of 
DM distribution. 

Since the observations at different off-GC angles are done dur- 
ing different time periods, to properly study the dependence of the 
line flux on the off-GC angle cj> one should take into account the 
time variability of the response of the SPI detectors. Several fac- 
tors have to be taken into account. First, the SPI instrument goes 
through a so-called "annealing" phase - heating of the detectors to 
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Figure 9. Dependence of the intensity of tiie positron anniliilation line at 
E = 511 keV on the off-GC angle. The solid line shows fit to the data in 
the fom const + iVe-'^^/'^'^^) . 



recover from a radiative damageF^Next. two of the 19 SPI detector 
have "died"13 The failed detectors also affect the response of their 
neighbors. To marginalize the effects of the changing response of 
the SPI detector, we split the entire data set into 7 periods, as shown 
on Fig.|7] The intervals are summarized in the Table|2] As both de- 
tector failures occurred soon after the end of an annealing phase, 
we chose to ignore revolutions 136 through 140 and 209 through 
215. The period 096-140 does not cover the essential part of the 
sky and therefore we skip it, leaving only 6 periods. 

For each of the periods, shown in the Table|2l we plot the dis- 
tribution of the line flux as a function of the off-GC angle (j). The 
results are summarized on Fig. [13] p. [16] One can see that none of 
these lines exhibits clear trend of decreasing from = 0° towards 
(j) = 180° . For each line (and each phase) we also compute the 
average flux F, standard deviation ap from the average, minimum 
(-Fmin) and maximum (Fmax). Our analysis shows that (a) 95— 100% 
of all points lie within Sap from the average (thus, the data is con- 
sistent with having flat spatial profile) and (b) the scatter of the 
data (Fmax — ^inin) IS much less than its mean value {F). Therefore, 
none of them cannot originate entirely from a DM decay. The cor- 
responding numbers for each line and each phase are summarized 
in TableU pagefTT] 

The positron annihilation line at i5 — 511 keV illustrates a 
situation, when a line of astrophysical origin is superimposed on 
top of the strong instrumental line. In this case, the data can be 
fitted by the constant, plus some function, depending on assumed 
shape of the source. Fig. |9] shows the dependence of the flux of the 
511 keV line on the off-GC angle of the SPI pointing. One can see 
that for the pointings with the off-GC angle less than 20° (about 
the size of the PCFOV of SPI) the 511 keV line flux contains a 
contribution from a sky source at the position of the GC, while for 
the pointing at larger off-GC angles the astrophysical source is not 
visible and the only contribution comes from the instrumental line, 
whose flux does not depend on the off-GC angle of the pointing. 



For details see SPI User Manual: 
|http://isdc.unige.ch/Instrument/spi/d oc/spi_um| 

Detector # 2 at revolution 140 and detector #17 at revolutions 214-215. 
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Figure 11. Life-time of the radiatively decaying DM as a function of the 
emitted photon energy. Region below the curve is excluded. 



5 RESULTS 

We analyzed the spectrum of SPI and found that none of the strong 
(i.e. detected with significance above 5a) lines can be interpreted 
as that of the decaying DM. This conclusion was based on the fact 
that variability of these lines over the sky is less than 10% (when 
moving from GC to the anti-center, see Fig[8](. At the same time for 
any realistic DM model such a variability would be greater than at 
least 60%. Thus, we exclude the possibility that one of the spectral 
lines, detected in the SPI background spectrum is a DM decay line. 

The non-detection of a DM decay line in the entire energy 
range of the SPI detector has enabled us to put an upper limit on the 
parameters of the DM particles. In particular, the Scr upper bound 
on the mixing angle of the sterile neutrino DM in the mass range 
40 keV - 7 MeV is shown on Fig.[TO] p.[T2] 

Our results are applicable to any decaying DM. To this end we 
also present the restrictions on the DM life-time (with respect to the 
radiative decay) as a function of the energy of emitted photon. The 
corresponding exclusion plot is shown on Fig.[TT] For example, the 
gravitino can decay into the neutrino and photon (similarly to the 
case of sterile neutrino) in supersymmetric theories with broken R- 
parity. Such an interaction is generate d via the loop effects (see e.g. 
iBorgani et aT]|l996l : iLola et al.ll2007l) . The restrictions on Fig.fTT] 
improve existing bounds on the lif e-time of such a gra vitino DM 
by several orders of magnitude (c.f. iBorganietaO 19961) . 

To present our results in the form less dependent on a particu- 
lar model of DM distribution in the MW, we show the 3a sensitivity 
towards the line search on Fig. [12] Note, that these results should 
be used with care, as the sensitivity depends on the assumed spatial 
profile of the source (because the effective area decreases with the 
off-axis angle, see discussion in Section|2]l. The results, presented 
on Fig.[T2]are valid for an extended source with the surface bright- 
ness which varies on the angular scales larger than (or comparable 
to) the size of the SPI fiel d of view (black solid line). T his plot is 
analogous to the Fig. 9 of iTeegarden & Watanab^ | |2006|) (TW06). 
However, a direct comparison of the Fig. 9 of TW06 and Fig.ll2lis 
not possible, since TW06 have assumed a different morphology of 
the extended source (10° Gaussian). Explicitly taking into account 
the dependence of the effective area of the SPI detector on the off- 
axis angle (see Section [2j, one can find that in order to make a 
direct comparison between the two figures, one has to "re-scale" 
the results of Fig. 9 of TW06 by an (energy dependent) factor of 
~ 1.5. This factor converts the sensitivity for the line, produced by 
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E [keV] 


Sign., (T 


AE [keV] 


Identification 


68.5 


11.4 


0.65 


66.7 Ge complex 


76.5 


58.8 


1.10 


75 Bi Ka 


87 


21.2 


0.90 


87 Bi Kp 


94 


5.2 


0.55 


91-105 GaZn 


134.5 


15.5 


0.90 


132-140 Ge complex 


143 


12.6 


1.20 


140-147 Ge complex 


177 


3.9 


0.95 


175 AsGe 


186.5 


7.3 


1.10 


184.6 GaZn 


193 


24.8 


0.75 


190-198 Ge complex 


200 


29.3 


0.60 


198-215 Ge complex 


205.5 


5.4 


0.50 


198-215 Ge complex 


240 


3.4 


1.10 


238 PbBi 


302 


3.9 


0.70 


301.5 GaZn 


311.5 


7.4 


1.10 


309.8 GaZn+K 


330.5 


3.2 


0.55 


328=?or331=PbTl 


385.5 


3.2 


1.45 


383=PbTl or ?? 


404.5 


5.5 


0.95 


403 Ga Zn-l-K 


F 431.5 


4.4 


0.55 




440.5 


12.4 


0.80 


438 ZnZn 


465 


4.3 


0.95 


470-485 NaNa 


1 51lt 


52.5 


1.25 


511 e+e" 


576 


5.2 


0.95 


574 GeGa 


585.5 


7.2 


1.10 


584.5 GeGa-HK 


597.5 


5.0 


1.45 


596-610 Ge complex 


754 


4.7 


0.85 


751 BiBi 


803.5 


3.7 


0.75 


803 BiPb 


812 


10.6 


0.95 


810CoFe 


819.5 


11.6 


0.85 


817CoFe+K 


827.5 


7.1 


0.75 


825 PbPb 


836 


12.1 


0.95 


834 MnCr 


845 


6.5 


1.30 


843 MgAl 


874 


8.8 


1.05 


872 GeGa 


884 


10.0 


1.05 


882 GeGa-t-K 


913 


5.2 


1.05 


911 AcTh 


937t 


3.1 


0.95 


?? or 935=MnCr 



E [keV] 


Sign., (T 


AE [keV] 


Identification 


990 


6.0 


0.85 


987 PbPb 


947t 


3.6 


1.50 




1014.5 


11.1 


1.40 


1014 MgAl 


1068.5 


3.3 


1.40 


?? or 1063=PbPb^ 


1079 


3.3 


1.00 


1077=GaZn 


1098 


6.8 


1.05 


1095=? 


1108.5 


19.9 


1.05 


1106 GeGa 


1118.5 


23.9 


1.05 


1117 GeGa-HK 


1127 


23.5 


0.85 


1124.5 ZnCu+K 


1234t 


4.5 


1.40 


1231 TaW 


1349.5 


3.5 


0.80 


1347 GeGa+K 


1368.5 


13.2 


1.55 


1368 NaMg 


1719.5t 


"3.3^ 


1.55 


1719 BiPb ^ 


1753.5 


4.0 


1.45 


?? or 1758=? 


1767.5 


11.1 


1.40 


1764 BiPb 


1781.5 


12.4 


1.45 


1778 AlSi 


1809+ 


15.2 


1.85 


1808 Mg 


1904 


3.2 


1.45 


1901=GeGa-HK or ?? 


2212 


6.7 


2.20 


2195-2223 BiPo, Al 


2225 


5.4 


1.10 


2223 HD 


2322 


3.1 


1.55 


2319=? 


2583.5 


3.5 


1.65 


2599=? or ?? ^ 


2616 


5.0 


1.90 


2614 PbTl 


2756 


9.5 


1.90 


2754 NaMg 


3002.5 


5.1 


2.45 


2993-3013 Al 


3176. 5t 


3.4 


1.95 


?? 


333lt 


3.4 


2.05 


?? 


3802 


3.6 


2.10" 


3800 GaZn-HK ^ 


4307.5 


3.6 


3.40 


4304 GaZn-HK 


4454 


4.2 


8.45 


4434 C 




3.5 


2.30 


?? ' 


5186.5+ 


3.5 


2.20 


?? 


5208. 5t 


4.0 


1.95 


?? 


5757+ 


3.4 


2.95 


?? 


6129 


20.5 


3.25 


6128.9 



Table 3. Lines, detected in the ON-OFF spectrum with the significance ^ Sct. Lines, marked with + cancel worse than by 90% in the ON-OFF spectrum 
(as compared with th eir flux in ON spe ctrum) and thus represent a "DM candidates ". The "Identification" c olumn indicates the probable identification of 
th e line in|W eidenspoi ntner et alj j2003l) . Lines, marked with "?" are not identified in lWeidenspointner et alj )2003l) . lines marked with "??" are not present 
in lWeidenspointner et alji2003l) . 



a source with the Gaussian surface brightness profile, into the one 
produced by a source of approximately constant surface brightness 
(see red curve in Fig. 1121). 



We have found a number of weak (with the significance 
3— 4(t) lines in the background-subtracted spectrum of SPI. These 
lines cancel by worse than 90% when subtracting OFF dataset 
(see Section |4l(. Apart from it we have found in the background- 
subtracted spectrum two lines with high significance - known lines 
at 511 keV and at 1809 keV. Any of these lines can in principle be 
a DM decay line. We analyzed each of them, by considering the 
profile of their intensity over the sky. Our analysis shows that none 
of these lines could be pure DM line (as their dependence on the 
off-GC angle does not show any clear trend to decrease towards the 
anti-center). The possibility that some of these lines are the super- 
position of instrumental and DM lines remains open. Quantitative 
analysis of the amount of DM flux admissible in a given line de- 
pends strongly on the model of the DM distribution in the Milky 
Way halo. Therefore it was not conducted here. 



6 DISCUSSION 

The purpose of this work was to understand how to search for the 
DM decay line with the SPI spectrometer and to check that none of 
the strong lines, present in the SPI background, was confused with 
the DM decay line. Our analysis shows that all the strong lines 
were, indeed, of instrumental origin and provides the upper bound 
on the flux of "weak" (3— 4a" above the background) lines, which 
leads to the corresponding restrictions (see Sec.[5](. To further im- 
prove the results, one needs to work with the weak lines (or lines, 
coinciding in position with instrumental ones). To do this one needs 
more sophisticated procedures of subtraction of the instrumental 
background (e.g. imaging). 

One of the most interesting cases of the coinciding instru- 
mental and celestial line is the positronium annihilation line at 
511 keV. An excess of positron annihilation emission on top of 
the strong instrumental line (related to positrons annihilating inside 
the detector) was noti ced long ago (for an incomplete set of refer- 
ences see e.g. Praiitzos 1 99 3; Milne etal. 1 999; Cheng etaL 1993; 
'Purcell et aU l 19971 ; iKnodlseder et alj|2005l ; IWeidenspointner et al.] 
12006 , 1 20071) . There exist many attempts of explanation of this ex- 
cess. In parti cular, it was attributed to the annihilating or decaying 
DM (see e.g. lBoehm et al.ll2004l ; [Hooper et alJl2004 lBoehm et all 
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Figure 10. Upper bound on the mixing angle of the DM sterile neutrino as a function of the sterile neutrino mass, obtained from the analysis of the background 
spectrum of the pointings toward the inner 13° of the Galaxy. For masses ^ 700 keV the restrictions from lYiiksel et all )2008l) (divided by the factor 2, due to 
the Ma jorana nature of the DM) are also shown in dashed line (see Discussion). For masses < 100 keV previous restrictions from HEAO-1 iBovarskv et alj 
l2006d) are also shown. The region above the curve is excluded. 
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Figure 12. 3(t sensitivity towards the line search in case of the flat dif- 
fuse source (thick black line). The results of lTeegarden & Watanabd j2006l 
Fig. 9), (rescaled to account for sensitivity towards the flat diffuse source, 
rather than 10° Gaussian) are shown in thin red line. 
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iFrere et alJboOTi : IPicciotto & Pospelovl l2005l : iRasera et al] 
The sterile neutrino DM with the mass rUs > 1 MeV pos- 



sesses decay channel Ns — > e^e v, with positrons annihilating 
either in flight or at rest, by forming the po sitronium atom (see e.g. 
iBeacom & Yuksell200^ : ISizun et alj200d) . Thus, it is possible that 
the decay of sterile neutrino DM contributes to such a line. The de- 
tailed analysis of this case will be reported separately. 

It should be also mentioned, that the region of masses between 
20 keV < ^ 40 keV remains inaccessible for the existing 

X-ray missions. The strongest restrictions in this region were pro - 
duced, using the data of HEAO-1 mission jBovarskv et al.ll2006cl) . 

When the work on this paper was at its final stage, the work 
of Y07 was published. Y07 obtained the restrictions on parameters 
of sterile neutrino in the range 40 keV - 700 keV. To facilitate the 
comparison, we plot the restrictions of ¥07 on Fig.[TO] (divided by 
the factor of 2 to translate them into the restrictions for the Majo- 
rana, rather than Dirac sterile neutrino DM, see footnote|7] p|6]l. As 
the data, used in our work, has about 5 times longer exposure than 
the INTEGRAL first years data, on which the results of Y07 are 
based, we could have expected results stronger by a factor ~ 2 in 
our case. However, the Fig. [TO] shows the opposite. The reason for 
this is as follows. For the SPI, the sensitivity towards the line search 
from a particular source depends on the shape of the source. In par- 
ticular, the results of TW06, on which the work of Y07 was based, 
were obtained under the assumption of a particular diffuse source 
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(10° Gaussian). As any realistic DM profile is much flatter than the 
10° Gaussian, the results of TW06 cannot be applied directly for 
the case of the DM line search. They should be rescaled to account 
for the diffuse nature of the DM source (c.f. Section|5}. Apart from 
this, the estimated DM signal from the inner part of the Galaxy is 
about 2 times stronger in Y07 than in our work. As the DM signal 
in the direction of the GC is the most uncertain, we have adopted 
the conservative flat profile everywhere inside the solar radius, to 
minimize this uncertainty. 
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Figure 13. Line flux as a function of the off-GC angle (p for tlie "candidate" lines from Table[3](page II U . For all lines the flux does not depend on the off-GC 
angle (with 95 — 100% of all points lying within 3 standard deviations from the average). Different colors represent different phases (see Table[2l p. [9). 
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Table 4. Characteristics of the spatial profiles of the candidate lines from Table[3] For each line (and for each of 6 phases) we compute the average (F), the 
standard deviation (average scatter of the points around its mean value) ap, minimal and maximal values and the ratio of (Fmax — ^miiO/^mim which gives 
the upper bound on the share of DM, present in the given line. 



